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ABSTRACT 

The production cross-section of the ^Pi charmonium state, can be predicted in 
Non-Relativistic QCD (NRQCD) using heavy-quark symmetry. We show that at the 
Large Hadron Collider a large cross-section for this resonance is predicted and it 
should be possible to look for the he through it decay into J/'ip + n even with the 
statistics that will be achieved within a few months of run-time at the LHC. 
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Non-Relativistic QCD (NRQCD) [Tj is an effective theory obtained from QCD useful 
for understanding the physics of quarkonia. In this effective description, states of 
momenta much larger than the heavy quark mass, m are excluded from the QCD La- 
grangian and new interaction terms are added to account for this exclusion. A cru- 
cial parameter is the relative velocity, v, of the quarks bound in a quarkonium state 
in terms of which the quarkonium state is expanded into Fock-components. It turns 
out that the QQ states appear in either colour-singlet or colour-octet configurations 
in this expansion where the colour-octet configuration evolves non-perturbatively 
into a physical colour-singlet state. The cross-section for the production of a quarko- 
nium H takes on the following factorised form: 

E -^(o?r+'i.)> (1) 

n={a,S,L,J} 

where F„'s are the short-distance coefficients, calculable in a pertubation theory in 
as, and On are local operators of naive dimension dn, describing the long-distance 
physics. The QQ pair produced in the short-distance process has a separation of 
a scale much smaller than 1/m which is potntlike on the scale of the quarkonium 
wavefunction, which is of order l/(mt>). The non-perturbative factor {O^) is pro- 
portional to the probability for a pointlike QQ pair in the state n to form a bound 
state H. The factorisation of the short-distance and long-distance parts of the cross- 
section guarantees the momentum-independence of the non-perturbative terms. 
These can be, therefore, obtained from one experiment at a given energy and used 
to compute the cross-section of the quarkonium state in a different experimental 
setting. 

Before this effective theory approach was developed, the production of quarko- 
nia was sought to be understood in terms of the colour-singlet model HZHS]]. While 
at lower energies this model was seen to provide an adequate description of the 
data, it was seen m in the phenomenology of large-p^ P-state charmonium produc- 
tion at the Tevatron [51 that colour-octet operators are very significant. Processes 
involving P-state quarkonia do not have a consistent description in terms of colour 
singlet operators alone [6J. Surprisingly, when data on direct J/^p production and 
on 'ip' production from the CDF experiment at the Tevatron was analysed, it was 
seen that it was necessary to include the colour-octet contributions for phenomeno- 
logical reasons [71, even though in the case of the S-waves the octet contributions 
are sub-leading in v. With the inclusion of the colour-octet contributions the full set 
of charmonium production data from the CDF could be described albeit at the in- 
clusion of non-calculable long-distance matrix elements 18,^. It was only the shape 
of the pT-distributions and not the absolute normalisations that was a prediction 
of NRQCD. Consequently, independent tests of NRQCD were necessary and sev- 
eral such proposals were made |[T0l[TTl[T2l[T3l[T4l[T5l[T6| . However, many of these 
proposals are not for large-pr quarkonium production and while they may be of 
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some phenomenological interest they do not provide a rigourous test of NRQCD 
because the NRQCD factorisation formula holds strictly only at large-pr- For a very 
comprehensive review of J / tjj production at the Tevatron and the related theory, see 
Ref. [17J. 

One interesting test of NRQCD comes from the study of the polarisation of J/ 'ip's 
at large-pr |T8|. The production of large-pr J/4''s proceeds primarily from the frag- 
mentation of single gluons and the QQ pair produced in the fragmentation process 
inherits the transverse polarisation of the gluon. The heavy-quark symmetry of 
NRQCD then comes into play in protecting this transverse polarisation in the non- 
perturbative evolution of the QQ pair into a J/ijj. The large-pr J/ip is, therefore, 
strongly transversely polarised. This is not true at even moderately low where 
the J/ ip is essentially unpolarised. The px dependence of the polarisation is, there- 
fore, a very good test of the theory ||T9| . 

The CDF experiment has measured the p^-dependence of the polarisation and 
they find no evidence for any transverse polarisation at large pt |201. Given the suc- 
cess of NRQCD in explaining the production cross-sections, this failure with respect 
to predicting the polarisation is, indeed, a shock. It may well be that the successful 
prediction of the production cross-sections of the various resonances was fortuitu- 
ous and that the effective theory is missing out on some aspect of the physics of 
quarkonium formation. It could be that the mass of the charm quark is not large 
enough to be treated in NRQCD. On the other hand, polarisation measurements are 
usually fraught with problems and it may well be that the problem is elsewhere. 
Finally, the problem may well have to do with the theoretical uncertainties in the 
prediction of polarisation. For example, the colour-singlet channel predicts the po- 
larisation of the J/ip to be longitudinal. So any effect that could substantially in- 
crease the colour-singlet contribution could change the full predictions of polarisa- 
tion quite drastically. To this end, a modified colour-singlet model with the produc- 
tion of J/ ^p's initiated by a scattering of a gluon with a Reggeized gluon has been 
considered [2T] but parts of the diffractive amplitudes involved in this calculation 
are not easily calculable. A more direct approach would be to study the effect of 
higher-order QCD corrections. These could substantially modify the theoretical ex- 
pectations regarding polarisation. Recent work 122| on NLO corrections to both the 
colour-singlet and colour-octet channels in the production of J/ ijj suggest that even 
these are not enough to understand the polarisation data. The situation is somewhat 
different in the case of T production [23J where the colour-singlet contribution, en- 
hanced by NLO and a part of the NNLO corrections, seems to be able to account for 
the data from Tevatron. For reviews of the current status of these calculations and 
their experimental consequences, see Ref. | |24ll25| . 

In this situation, it is worthwhile looking for other tests of NRQCD which suc- 
cessfully navigate between low-px and polarisation. Such a suggestion had been 
made years ago in the context of charmonium production at Tevatron 1126 II : the 
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production of he, the ^Pi charmonium state. In NRQCD, this state is produced in 
the colour-singlet mode and through the production of an intermediate octet ^Sq 
state. The non-perturbative matrix element for the transition of this octet state to 
the physical ^Pi state can be inferred from other non-perturbative parameters fixed 
at the Tevatron. This is a consequence of the heavy-quark symmetry of NRQCD. 
Consequently, one can predict the rate for he production in NRQCD. In this letter, 
we investigate this prediction in the context of the Large Hadron Collider (LHC). 
One channel which may be suitable for the detection of the he is its decay into a 
J/il) + Tx. The decay branching fraction for he J/V^ + vr has been estimated from 
spectroscopy. 

It may be argued that the measurement of the other charmonium resonances like 
the J /if), x's and the if)' will already provide the tests of the NRQCD factorisation 
formula. The non-perturbative parameters have been determined at the Tevatron 
and the factorisation formula implies that these are not momentum-dependent. So 
it should be possible to predict the cross-sections for these resonances at the LHC 
and check for the validity of NRQCD. While this is true, it must be remembered 
that for several years quarkonium production has also been studied in terms of 
a phenomenological model known as the semi-local duality model or the colour- 
evaporation model [27J. In this model, it is assumed that the open-charm cross- 
section integrated over the region between 2m and the open charm threshold should 
be equal to the sum of the resonance cross-sections. The resonance cross-section is 
then some fraction of the open charm cross-section integrated over this mass range. 
The fraction is unknown apriori but is fixed by comparing to the data - it is the 
analog of the non-perturbative parameter that appears in NRQCD computations of 
the cross-section. This approach is seen to provide a reasonable description of the 
data from Tevatron 11281 129II . However, it must be borne in mind that the separa- 
tion into perturbative and non-perturbative parts in this model is not rigourously 
provided by a factorisation formula as in NRQCD and, consequently, the fractions 
(non-perturbative parameters) that are determined by fitting to Tevatron data are 
not guaranteed to be energy-independent. If the energy dependence of these pa- 
rameters is large, then it will not be possible to use the semi-local duality approach 
in any predictive way at the LHC. However, it may so happen that, in actual prac- 
tice, the energy dependence of the fractions turns out to be small in which case 
semi-local duality will be able to predict the resonance cross-sections at the LHC 
as well as NRQCD can. But these predictions, in the semi-local duality approach 
can be made for only those resonances which have been measured at the Tevatron. 
It is not possible to predict the cross-section for particles which have not been de- 
tected at the Tevatron within this model approach. The search for he at the LHC is, 
therefore, important in establishing NRQCD as the correct theory of quarkonium 
production. 

It is also worth emphasising that the he had eluded experiments for a long time 



3 



and it is only recently that its existence has been verified in e+e experiments at 
the CLEO Il30ll . The LHC is expected to produce this resonance copiously and it 
may provide a study of this resonance in various decay channels and may help 
understand its properties. 

At the LHC, the production of the he proceeds through the following partonic 
subprocesses: 

g + g 

g + g ^'S^^^+g, 

q + q ^'St^+g. (2) 
The large-pT hadronic production cross-section is given as 

dpT 

^ f dy f dxiXiGa/pixi)x2Gb/pix2)- — 3^ — y^^^^ ~^ "^^"^^Ljc). (3) 

In the above expression, the sum runs over all the initial partons contributing to the 
subprocesses; Ga/p and Gb/p are the distributions of partons a and h in the hadrons 
with momentum fractions xi and X2, respectively The expressions for the singlet 



and the octet subprocess cross-sections, da/dt, are given in Refs. ||3111 and BH], re- 
spectively 

The ^Sq'^ he is mediated by a gluon emission in a El transition. To fully 
determine the production rate we need the colour-singlet matrix element for the 
^Pi state (^Oi''{^Pi)'j and the value for the colour-octet matrix element that takes the 

octet ^5*0 state to a h^, (|(98^ (^S'o)^. The colour-singlet matrix element is related to the 
derivative of the wavefunction of at the origin by 

otCPi)) = ^\R'm'- (4) 

The Tevatron data on Xc production fixes [9] the colour-octet matrix element which 
specifies the transition of a ^6*1 octet state into a ^Pj state. We would expect from 
heavy-quark spin symmetry of the NRQCD Lagrangian that the matrix-element for 

^5*0^^ he should be of the same order as that for ^sf^ — > ^Pi. This is because the 
essential difference between these transitions comes through the magnetic quantum 
number so that the corrections to this equality will be of 0{v'^) ~ 30%. For the 
derivative of the wave-function we use a similar argument to fix it to be the same as 
for the Xc states. 
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Figure 1: The cross-section for he production as a function ofpx cut for different choices of 
QCD scale 

With these inputs, we have computed the cross-section for he production in pp 
collisions at the LHC {^/s = 14 TeV). We have computed the cross-section integrated 
over pt with a lower p^ cut. In Fig. 1, we present the results for the pr-integrated 
cross-section as a funtion of the pr-cut for three different choices of the QCD scale: 
Q = Mr/2, Mt and 2Mt. We have used the CTEQ 4M parton densities [32] . The 
cross-section has been folded in with the branching ratio of the ^Pi state into J/tp + n 
and the J/ip l^l^, where / = e or /i. We have integrated over the rapidity interval 
—2.0 < y < 2.0. For the singlet matrix element, we use the value extracted from Xc 
decays, which is (^Oi''{^Pi)'^ = 0.32 ||33| and for the octet matrix element we have 

(C»8^^(^5o)) = 0.0098 0. With these inputs, we find that the cross-section for h^ 
production (folded in with the decay fraction into a J/tjj and vr, which we take to be 
0.5% Il34l and a 6% leptonic decay branching fraction of the J/ip) is large enough to 
have a substantial number of events with the statistics that will be acquired in the 
first few months of LHC running. Varying the QCD scale between the largest and 
the smallest values that it can take results in a variation in the cross-section which 
is about a factor of 2. While the results for the cross-section for he production in 
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Figure 2: The cross-section for he production as a function ofpT cut for the range of allowed 
values of the octet matrix element 



Fig. 1 show the variation with respect to QCD scale inputs, in Fig.2 we display the 
uncertainty in the cross-section coming from varying the value of the octet matrix 
element. We expect a 30% variation about the central value of 0.0098 for the octet 
matrix element. The two curves in Fig. 2 correspond to the upper and lower values 
that the octet matrix element can take. In Fig. 2 the QCD scale is taken to be Mp. 
The variation in the cross-section due to the change in the octet matrix element is 
about 60%. 

In Fig. 3, we show the pr-integrated cross-section choosing different parton 
density sets. In addition to the CTEQ 4M densities used earlier, we use the LO 
CTEQ 1132 II and GRV densities ||35| . It is only at low values of px that a sizeable 
change in the cross-section due to the variation of the parton density inputs can be 
seen and even at a pt value of 10 GeV the variation is not more than about 30%. 
The decay branching fraction of he into a J/tp + n could be as large as 1% |34|, and 
if we use this instead of the 0.5% used in the above calcuations we could have a 
production cross-section which is twice as large. 

The Pt distribution Bda / dpx is shown in Fig. 4. We have plotted the octet and 
the singlet contributions separately. We find that, over a whole range of large px, 



6 




0.001 I 1 ' ' ' 1 1 

20 40 60 80 100 120 

p_T cut 

Figure 3: The cross-section for he production as a function of px cut for different parton 
distribution sets 

the singlet contribution is negligible and that the he is produced almost exclusively 
from the colour-octet channel. 

We would like to conclude by making the following points: 

• In spite of the success that we have had in understanding charmonium pro- 
duction at the Tevatron using NRQCD, we still need to have independent 
tests of this effective theory because the colour-octet parameters, and conse- 
quently the normalisations of the cross-sections of the various charmonium 
resonances, are not given by the theory but only fixed by fitting to the data. 

• Polarisation predictions for J/ ^ and tjj' at large=pT/ considered to be good tests 
of NRQCD, disagree violently with what is measured by the CDF experiment 
at the Tevatron. 

• The production of he in NRQCD is a good test of the theory because: i) it is 
a prediction for large-pT production where NRQCD factorisation is expected 
to hold, ii) the cross-section can be predicted because the relevant colour-octet 
parameter can be inferred from octet parameters measured in x production 
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Figure 4: The pt distributions for production 

at the Tevatron and using heavy-quark symmetry and, iii) the cross-section 
is very large at the LHC and should lead to easy detection of the resonance. 
Moreover, the cross-section measurement is much simpler than measuring the 
polarisation of the charmonium state. 

• Such a prediction for the cross-section of he production cannot be made in the 
alternative approach to quarkonium production, viz., the semi-local duality 
model 1271 1^ in . If the prediction of NRQCD is verified it will certainly 
establish it as the correct approach to quarkonium physics. 

In conclusion, even with the statistics accumulated with a few months of LHC 
running the charmonium resonance, h^., can not only be detected but its properties 
can be studied in detail. We have presented predictions of NRQCD for the pro- 
duction cross-section of the he and so the study of this state at the LHC will help 
test NRQCD independently and provide us more understanding of the physics of 
quarkonium formation. 
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